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Abstract 

An interesting situation occurs when the linearized dynamics of the 
shape of a formally stable Hamiltonian relative equilibrium at nongeneric 
momentum 1:1 resonates with a frequency of the relative equilibrium's 
generator. In this case some of the shape variables couple to the group 
variables to first order in the momentum perturbation, and the first order 
perturbation theory implies that the relative equilibrium slowly changes 
orientation in the same way that a charged particle with magnetic moment 
moves on a sphere under the infiuence of a radial magnetic monopole. In 
the course of showing this a normal form is constructed for linearizations 
of relative equilibria and for Hamiltonians near group orbits of relative 
equilibria. 
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1 Introduction 



Suppose we are given a symplectic manifold (P, lu) and a compact Lie group G 
acting freely and symplectically on P with Ad*-equivariant momentum mapping 
J . P ^ 2*. Let the Lie algebra of G be g and H : P R he a. G-invariant 
Hamiltonian. Let us fix our attention on a specific relative equilibrium p^, with 
generator ^e, so that the integral curve of the Hamiltonian vector field Xh at 
Pe is e'xp{£,et)pe- Let J{pe) = Me, let Gfj_^ be the isotropy group of fie under 
the coadjoint action of G on g*, and let g^^ be the Lie algebra of G^^. For 



more information on these basic definitions, see, for example, Abraham ancj 
MarsdcI^ |197S| ], |Marsdcn| [ |1992[ , or [Marsdcn and Ratirj |1994[ . 

Suppose that G = 50(3), although the results of this article are much more 
general than that. Since the action is free, a neighborhood of the group orbit 
Gpe may be (with respect to the left action) equivariantly projected to SO(3) 
itself. Let us realize S0{3) as the unit circle bundle of the 2-sphere in M'^, ar- 
ranging that the orbit exp(^t)pe be projected to the the unit circle in the tangent 
space of the vertical vector k. The system, when at the relative equilibrium pe, 
and when viewed through this projection, appears merely as a point, say V, ro- 
tating uniformly on the unit circle in the tangent space of S'^ at k, and the other 
relative equilibria gpe, g € 150(3), appear as this same motion but reoriented 
by the rotation g. Suppose Pe is formally stable, and suppose the system is then 
perturbed from the relative equilibrium Pf,. If /ie 7^ then the motion of the 
system is bound to remain near exp(^eR)pe, as shown in Patrick [ 1992[ (exten- 
sions to the nonfree case may be found in prtega and Ratiti 1 1999| ] and Lerman 



and Singer | 1998( |). The perturbation affects the motion of the point V by su 



perimposing on its original motion some small vibration; however the motion of 
V does not ever carry that point far from its original circular path. 

The situation is different when /ig = 0, for then Gf^^ — 50(3), and under 
perturbation the system is bound only to remain near Gpe- The effect is that, 
after a perturbation from pe, the point V still moves nearly circularly in tangent 
spaces of 5^, but the base point of the tangent space slowly changes location on 
5^. In fact, as shown in Patrick [1995|, the base point B moves as a point charge 
under the influence of a radial magnetic monopole. This result is obtained by 
constructing coordinates which make an open neighborhood of Gpe in P into 
an open neighborhood of Z(r*S'0(3)) x {p^} in T*5'0(3) x P^^, where P^^ = 
J~^(/Xe)/G^^ is the symplectic reduced space, pe € P^, is the equilibrium of the 
reduced space corresponding to Pe , and Z(T*SO{3)) is the zero section of r*S'0. 
In the new coordinates, to first order in the momentum perturbation, the "group 
variables" in r*50(3) and the "reduced variables" in P^^ become decoupled; 
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truncating the higher order interaction terms gives a new Hamiltonian system 
on T*SO{3) X P^^ called the drift system. The drift system has an additional 
normal form symmetry. The Hamiltonian system for the slow motion of the base 
point B is obtained from the drift system by ignoring the variables on P^^ and 
reducing the resulting Hamiltonian on T*SO{S) by the normal form symmetry. 
The approximating, truncated Hamiltonian of the drift system may be obtained 
merely by calculating the nilpotent part of the linearization of the vector field 
Xh/:^ (p) = Xh ip) — £,eP at Pe , where ^p = (^) p (p) denotes the infinitesimal 
generator of at pe . 

Actually, the results in Patrick | 1995 | depend on a nonresonance condition 
that the frequency |^e| cannot equal any linearized frequency of the reduced 
system at pe- This article considers the problem of constructing a Hamiltonian 
system for the motion of the base point B in the case where there is such a 1:1 
group-reduced resonance. The main result of this article is that, in the presence 
of resonance, some of the reduced and group variables couple to first order 
in the momentum perturbation, with the result that B moves as a charged 
particle with magnetic moment on a sphere under the influence of a radial 
magnetic monopole. However, this "particle" regularly exchanges its charge 
with its magnetic moment. 

The construction of the Hamiltonian system which models the motion of 
B in the 1:1 resonant case generalizes and parallels the nonresonant case. In 
the course of the construction the linearization of Xh^ at Pe is put into a 
certain block normal form through a Witt or Moncreif decomposition of the 
tangent space Tp^P. The objective of this normal form is the separation of 
the group and reduced motions and it is not the full normal form of the lin- 
earization as a infinitesimal symplectic operator. The linear normal form is 
then extended to a neighborhood of the group orbit Gpe using the the isotropic 
embedding/equivariant Darboux theorem. These normal forms — for lineariza- 
tions of relative equilibria and for Hamiltonian systems near group o rbits of 
relative equilibria — are the second main contribution of this article. In Patrick 



and Roberts 1999 these normal forms have been used to study the structure 



of the set of relative equilibria. 

A new aspect that is absent in the nonresonant case but emerges in the 
resonant case is the appearance of a gauge group. There is too much freedom 
inherent in the normal forms to exactly fit the needs of first order agreement 
between the original and drift systems. This slack is taken up by the group 
(R, -t-) and presents itself as an inherent freedom in the choice of the normal 
forms. This gauge freedom can be used to simplify the drift system, a motif 
which is quite useful in this work. 

The Hamiltonian system modeling the motion of B is approximate since 
it is the result of a truncation, so there arises the question of how well this 
model Hamiltonian reflects the behavior of the real one. To check this I have 
numerically simulated the system of two (identical) axially symmetric rods 
which are joined by a frictionless ball-and-socket joint, but otherwise move 
freely (Patrick |1989, 1991 1). The relative equilibria of this system correspond 
to motions such that the two rods spin on their axes while otherwise main- 
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taining constant mutual orientations, while the whole assemblage rotates about 
some fixed axis. There exist, due to the system's multiple rotating parts, rela- 
tive equilibria which are formally stable and which have zero (hence nongeneric) 
total-angular-momentum. In fact there is a continuum of relative equilibria with 
zero total-angular-momentum, and it happens that parameters may be chosen 
so that there are relative equilibria at zero total-angular-momentum and at the 
1:1 group-reduced resonance. The third main contribution of this article is the 
actual calculation of the drift system for the coupled rod system near one such 
relative equilibria, and the verification of the drift approximation by comparison 
of simulations of the coupled rod system itself and certain predictions of its drift 
system. 

Here is an overview of this work. I begin in Section 2.1 with an analysis 
of the linearization dXn^^ (pe) of Xh^^ at Pe, by splitting the linearization into 
semisimple and nilpotent parts: dXH^^ipe) = Sp^ + Np^. The key tool in the 
analysis is a certain Moncrief or Witt decomposition of Tp^P: 



Whereas this decomposition is obtained in Patrick | 1995| ] for the nonresonant 



case with the aid of the generalized eigenspaces of dXH^^{Pe)^ and so was 
dXn^ (pe)-invariant, here the decomposition is only S'p^ -invariant. Using this 
decomposition, the nilpotent part Np^ acquires a block form which gives rise to 
certain operators: 

In the nonresonant case iV^J = and N^^ = necessarily. I expose, in Sec- 
tion 2.2, the special properties of these operators, the most notable of which are 
that Np'^ is symmetric, N^'^ and N^^ are dual, and that these operators have 
certain commutation relations with operators adj^ and coad^^ . 

Next, in Section 3, I focus on the case where the reduced and group spec- 
tra have intersection {±iAp^}, a purely imaginary eigenvalue and its conjugate, 
each of which occur with multiplicity 1 in both the reduced and group spec- 
trum. The generalized eigenspace of the reduced linearization corresponding 
to {zbiAp^} is symplectic and has dimension 2, and so is linearly symplecto- 
morphic to (K^ = {xi, 2:2}, dxi A dx2), and the operator N^^ is zero on the 
sum of the complimentary generalized eigenspaces. Consequently, the splitting 
of Tp^P refines and the operator N^^ may be replaced by another operator 
Np^"^ : Qfj^^. After this simplification, the splitting of Tp^P and the 



isotropic embedding theorem (Marsden |1981],Weinstcin |1977|) together give 
a map which transforms the Hamiltonian H, to first order near G^^^pe, to the 
Hamiltonian 

near {0} x Z{T*Gf,J in the phase space x r*G^^, where Z{T*G^J denotes 
the zero section. In addition to the expected invariance under the left action of 
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G^^ , Hdrih is invariant under a diagonal action of the toral subgroup generated 
by ^e! this is the normal form symmetry. 

In Section 4 I further assume that G is the largest compact continuous sym- 
metry group for an ordinary mechanical system, namely 5*0(3) x (S*^)". After 
various manipulations, including an Abelian reduction, the drift Hamiltonian 
becomes 

-ffdrift = (tTi^ + TT2^) + ^hT^z' + k{'KiXi + 7120:2) + aTTs 

on the phase space x T*SO{2i) = {{x, A, tt)}. Here Ii and I2 come from N"^"^, 
K, a single coupling constant, is all that remains of N^^^, a is a constant, and 
there are two possibilities: the relative equilibrium can be either a "+" type 
or a "— " type. In this context the gauge freedom has the effect of making Ii 
arbitrary and interest is focused on the dynamics near to a; = 0, tt = 0. The 
Hamiltonian ifdrift is defined on a phase space of dimension 8 and has symmetry 
5*0(3) X S^, and so defines a system which is completely integrable. Although 
a complete analysis in the general case seems difficult, I analyze this system 
in Section 4.1, where I show that the system has, for example, some singular 
reduced phase spaces and spectrally unstable relative equilibria with homoclinic 
connections. It is in Section 4.2 that I show that the drift system reduced by the 
normal form symmetry can be cast as a charged particle with magnetic moment 
moving on the sphere while under the influence of a magnetic monopole. 

Finally, in Section 5, I numerically investigate the dynamics of the coupled 
rod system near one particular resonant relative equilibrium and compare this 
dynamics with that of the drift system. I explore three distinct regions of phase 
space: 1) within zero total-angular-momentum, 2) near a stable relative equi- 
libria of the drift system, and 3) near a spectrally unstable relative equilibrium 
at a singularity of the drift system. The comparison of the two systems is ham- 
pered by the implicit nature of the coordinates relating them, but in the first 
two comparisons agreement between the two systems is obtained uneventfully. 
However, the singular points of the third comparison involve an unexpected 
reconstuction phase jump; the situation is delicate and small perturbations are 
required to elicit quantitative agreement between the two systems. Nevertheless 
it becomes clear that many elements of the dynamics of the coupled rod system 
near the resonant relative equilibrium are indeed captured by the drift system. 

2 The linearization 

Here are the basic notations: 

1. Pp is a relative equilibrium, is the generator of and the momentum 

of Pe is /ie = J{pe). 

2. G^^ is the isotropy group of He under the coadjoint action of G on g*, 
and fly^e is the Lie algebra of G^^. CoAdg = (Ad^-i)*, g G G, and 
coadj = -(adj)*, ^ e g. 
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3. H^^ = H — J^^, so that H^^ has a critical point at Pe- The Hessian of i?^^ 
at Pe will be denoted by (PH^^ (jpe) and the linearization of ^i/^^ at Pe will 
be denoted dXn^^ (Pe)- 

4. Without loss of generality, since everything is local near to Pe and pe is 
regular, the Marsden-Weinstein symplectic reduced space (P^^,w^J ex- 
ists; let TT^^ : J~^{fie) Pfj.^ be the projection. The reduced system has 
this phase space with Hamiltonian i/^^ defined by H^^ir^^ = H\J~^{fj,e). 
Also, Pe = 7r^_.(pe), and the linearization of the reduced system at pe is 
denoted by cLXh^,^ (Pe)- Suppose that Pe is regular, which means S.Pe 7^ 
for all C G 0- 

5. The reduced spectrum is the spectrum of the linearization dXn (pe) at 
Pe of the reduced system. The group spectrum is the the spectrum of 
adj. : ^ fl). 

Let dXH/:^ (pe) = Np^ + Sp^ be the Jordan decomposition of dX^i^^ (Pe) into 
its semisimple part Sp^ and nilpotent part Np^ . The aim of this section is an 
analysis of dXn^^ (Pe), focusing on its nilpotent part Np^ . Do not assume that 
the reduced spectrum and the group spectrum are disjoint. 



2.1 The splitting of Tp^P 



I begin by deriving a "Moncreif" splitting (Marsden |1981]) of Tp^P which is 
slightly weaker than its analogue in Patrick |1995|, in that the "reduced" part 
(below M^ed) of the splitting is not (and cannot, in general) be constructed to 
be dX}j^^ (Pe) invariant. The details are similar to those in Patrick | 1995{ ; I will 
not belabor them here. The subspace kerdJ(pe) is -invariant, so one can 
choose an Sp^. -invariant subspace W^cd such that 

ker dJ {pe) = T^od ® fl^ePe- 



The subspace W^od is symplectic, and T-K^^\Wred ■ M^ed — > Tp^ 
symplectomorphism. Choose an Ad-invariant complement b to Qn^ , so that 



is a linear 



= b ® 0^^. 

The subspace bpe is symplectic, and C {Wrcd ® bpe)"""" as an S'p^^ -invariant 
Lagrangian subspace. Choose an S'p^ -invariant Lagrangian complement^ Z to 
fl/je in {Wred ® faPe)""^, giving the S'p^ -invariant splitting 



Tp^P = W^rcd © St^^Pe ®Z®bpe 



(1) 



As already stated, Ttt^^ is a symplectomorphism between W^ed and Tp^P^^. 
Also, dJ{pe) is a linear isomorphism between Z and g*^, and a symplecto- 
morphism between bpe and the tangent space T^^(G^e) at /ig to the coadjoint 

^ Since Sp^ is semisimple, every Sp^ -invariant Lagrangian subspace has an 5p^ -invariant 
Lagrangian complement. This is a general f act: the proof is a simple modification of the 



argument found at the bottom of page 401 of Patricls [1995 
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orbit Gfie- Thus becomes, through these identifications, the iSp^ -invariant 
sphtting. 



T- P 



(2) 



With respect to the decomposition (||), let dXH^^iVe), Sp^, and Np^ have 
blocks [A^^], [S**^], and [iV*-^], respectively; the block form of aj(pe) becomes 












-Id 







Id 












(3) 



where ujp,^ is the reduced symplectic form of P^^ and uj^^ is the Kostant-Souriau 
form on G/Ltg. In (|]), the subspaces g^^Pe, &Pe, and ker (iJ(pe) = W^od ffi Qti^Pe 
are dXjj^^ (j)e)-invariant, and hence are A'p^ -invariant, and this implies certai n 
of the A^^ and iVp-^ vanish. Using the identities (Proposition 5 of Patrick 1 1995f ) 



Xh^^ (Pe) = - coad^^ d.J (pe) 
dXn^^ {Pe)v{Pe) = -(ad^^ v)p{Pe), V & 9 

to calculate the diagonal blocks of dXH^^{Pe), the Jordan decomposition be- 
comes 



dXH!:APe) = 



d^H^^Pe) 








ad^ 







^13 

ad^^ 



5" 














g22 




























544 








ad 








iVl3 





Pe 


^22 

Pe 


^23 

Pe 











^Pe' 














^44 

Pe 



Since Sp^ and Np^ commute, so do and 
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and S*^^ is semisimple and N^^ 
is nilpotent since Sp^ and Np^ are. Thus, dXu^^ (pe) = 5" + is the Jordan 
decomposition of of the reduced linearization dXH^^{pe), which is semisimple, 
since Pe is formally stable. Consequently, — dXn^^ (pe) and N^^ = 0. 
Similarly, N^^ = N^^ = N^j = 0, and 



adf 



533 = adj 



S^^ = ad? 



Thus, 



d^H^SpA) 








— ad^^ 

















(ad^J- 



(4) 
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and 



N, 
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N^l iV2f 



(5) 



As will be seen immediately below, a large amount of information can be dis- 
cerned from structure of Np^ and Sp^ visible in (^) and (||) . 



2.2 Properties of Np^ and its subblocks 

Directly from (|^), (Np^)^ = so that Np^ has nilpotent order at most 3. This 
upper bound on the nilpotent order occurs because dX}{i^^{Pe) is semisimple 
when suitably restricted and projected to the summands of (|l|), so that Np^ 
must "transport" between those summands. Moreover, the "direction" of this 
transport is "one way", since ker c?J(pe) and Q^^Pe are invariant subspaces, and 
the "transport" stops at Qfi^Pe, on which Np^ is zero: 

Image iVp^ C kerdJ(pe), (6) 
Np^ (ker dJ (pe)) ^ Qt^^Pe, 

Npj9^^Pe = 0. (7) 

When iteratively acted upon by Np^, a vector can make just 2 stops before 
annihilation: ker c?J(pe) and g^^Pe! hence the nilpotent order of Np^ is a most 3. 
By (0), 

kerTp^TT^^ = g^^Pe C ker(A'pJ kerdJ(pe)) , 
so : Tp^Pp^^ ^ is the unique linear map such that 

kerdJ(pe) > 9fj^^Pe ^^"^^j QtJ.^ 

Tp^Pfj,^ 



(8) 



Since this diagram does not depend on choices made in the construction of the 
splitting (|l|), neither does N"^"^ . Also, N'^^ may be calculated merely by reference 
to the diagram; it is not necessary to calculate all summands of (^). Similarly, 
by (|), and since, by (||), 

kerdJ(pe) C ker(r7r^^ A'pJ, 
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: 0*^ — + Tp^P^^ is the unique Hnear map such that 




T- P 

Since dXui^^iPe) is infinitesimaUy symplectic, so is Np^, so that 

{Npyuj{p,)+u{p,)Np^=0, (9) 

and inserting (||) and (||) into (^) gives 

Nl] = -{Niycol^, [Niy^Nll (10) 

so A'^^ is symmetric, and in a certain sense, Np'^ and Np^ are dual. Similarly, 
using and (j^) to write out what it means for Sp^ and Np^ to commute, 

dXn,^ {Pe)Nll = -Nil coad^^, (11) 
NlldXH^^ (Pe) = - ad^^ TVp^i , (12) 
ad^^ Nil = Ar23 ^^^^^^ (^3) 

These properties will yield the important normal form symmetry of the drift 
system. 

Temporarily let pr2 be the projection onto the second factor of ker (iJ(pe) = 
W^cd©0M=Pe- Then 

ker(dJ(pe)|(M4cd + bpeT^) - Q^.,Pe C ker(pr2 Np^{W,,d + bpe)"^), 

so : Qi^^ X Qi^^ — > M is the unique bilinear from such that 
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Thus, iV^^ may be calculated without calculating the Lagrangian complement 
Z of (§. 

However, A'^^ does depend on the particular choice or W^ed, even if W{-ed 
is an S'p^ -invariant complement of Sfj,ePe ker (iJ(pe)- To see this dependence, 
choose a linear map A : H^od — > 9^,^ , which commutes with Sp^ , and then choose 
a new Wred where 

W^cd ^{w + iAw)pe \weW,,d}- (15) 

Also, choose a map B : b which intertwines the adjoint action and then 

choose a new b' by 

b' ^ {ri + Bri\7j eb}. (16) 

Given v G g*^, the process of calculating {Npf^v, given by the analog of ( p^ 
for the new choices, is 

1. pick z' e {W^^^ + b'pe)""^ such that dJ{pe)z' = v] 

2. calculate Np^z'; 

3. project Np^z' to Qfi^Pe = S/j^ using the splitting of ker (iJ(pe) defined by 

with result {N^ffi^. 

It is straightforward from (^), (p^), and (p^), that 
so given v £ 0* ^ , let 



z' = © © ® ujl^B*iy. 



Then 



iVj^^z' = iVpif © {N^lul^A*v + {Niyv) © © 



13.,^ AAjl\ 

21, i /I*,, /lArl3,, , /Ar23\b, 



iVp^> © AA^ji^ © © 



+ © (Af;>»,^ A*z/ - AN^y + (iv;f ) V) © o © o, 

which by ( p^ is the appropriate decomposition of Np^z' . Thus, using (10), 

{Nlf fv = N'pHA*. - AN'^y + 

= -{{N;,yA*+AN;,^),y+{N^y^. 

This freedom to adjust N^^ by adjusting W^cd will give the important gauge 
freedom of the drift system. 

This gauge freedom can be encoded as a group, as follows. Temporarily set 
Epe = Tp^Pfj^^ X 0^^^ X 0*_^ X Tf^^Gp^ , and fix one particular splitting of type (^, 
giving a symplectomorphism, say 0o : Tp^ P — > Ep^ . Given A : M^cd — 0/^^ and 
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B : b ^ 2fi^ as above, there is another symplectomorphism (j)A.B ■ Tp^P — > E; 
and hence a unique Aa,b in the symplectic group Sp{Ep^ ) such that 




Ep. > Ep. 

Conversely, A is determined by Aa,b, since if a; G Tp^P^^ , ^ e g^^., and w G M^cd 
is such that (l)oiw) — x, then 

AA,bix) = (t>A,B{w) = (f)A.B{w + {Aw)pe - {Aw)pe) 

= TlTf,^ (w + {Aw)pe) ® {-Aw) =x® {-Aw). 

Similarly B is determined by Aa,b, so the pairs {A,B) have a natural group 
structure given by the injection {A, B) i-^ (f'A,B, and a simple calculation gives 

4>A+A\B+B' = 4'A,B ° 't'A'.B'- 

Thus, at this level, the gauge group is the additive Abelian group of pairs {A, B) 
such that A commutes with Sp^ and B commutes with ad^^ . 

3 The simplest resonant case 

All the above has been established without any particular presumption on the 
way that the group and reduced spectra intersect. There are some apriori fea- 
tures of the spectra: the spectra are purely imaginary (or zero) and are invariant 
under change of sign. Since pe is presumed to be formally stable, zero cannot 
occur in the reduced spectrum, and hence not in the intersection of the spectra. 
Thus, the following resonance condition is the simplest possible case beyond an 
empty intersection: 

The intersection of the spectrum of ad^^ and the spectrum of 
dX}j^^{Pf.) is {zbiAp^}, \p^ > 0, and each of the eigenvalues 
zbiAp^ occur in each of these operators with multiplicity one. 

While many of the results below remain true or have analogues in a more relaxed 
environment, for simplicity, this resonance condition will be assumed for the 
remainder of this article. 

3.1 The linearization 

Given the resonance assumption, it is natural to consider certain spectral split- 
tings of the factors Wrcd and g^^ in In particular, let W^Jj be the generalized 
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eigenspace of ±iXp^ for dXn^^ {Pe), and W^^^ be the sum of the other general- 
ized eigenspaces, so that W^J^ is two dimensional, symplectic, and there is the 
symplectic splitting 



(18) 



Similarly, let be the generalized eigenspace of ±iAp^ for the linear map 
adj^, and gj^^ be the sum of the other generahzed eigenspaces, so g^^ is two 
dimensional, and 



=0L®0L- (19) 

Of course, there is also then the dual splitting g* ^ = g°* ® g^* . 

So we have a refinement of the splitting (H) into "nonresonant" (with super- 
script 0) and "resonant" (superscript 1, dimension 2) parts: 



-,0 



'0 



(20) 



As one might expect, N^'^ and A'^j^ localize to the resonant parts: 

C ker Nil , Image C gi ^ , 
g^l^ C ker N^^ , Image iVpif C W,l^ . 



(21) 
(22) 



Indeed, to show the first of (^, it suffices to find a subspace, say 14^°^ of 
ker(iJ(pe) such that 



-M^.;id and N,^W^^^=Q. 



(23) 



The result then follows by using to reverse the vertical arrow of (^) . The 
first of (23) can be assured by setting W^^j to be the sum of the of nonres- 
onant (i.e. not zbiAp^) generalized eigenspaces of dXn^^iPe), and the second 
of ( p3| ) follows since is dA//^^ (pe) -invariant and that dXHi:^{jpe)\W^^^ is 

semisimple. For the second of (^T]), Np^ maps the ±iAp^ generalized eigenspace 
of dXff^^ (pe) to itself, since dXn^^ (pe) does that. Thus the image of Np^ is 
contained in the intersection of that generalized eigenspace with Q^^Pe, which is 
exactly Q^^Pe- From ( |2l| ) one gets ( p2[ ) by using the duality (|l^ between A^i 
and iVpi3. ' 

With respect to (pO|), and in view of (|), (1^) and (H), the hnear map iVp^ 
has the form 



Np, = 



























^131 

Pe 








^211 







^23 

Pe 
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where N^^^ = A^pJ|W^ed ^^'^ where A''^^^^ is iV^^^ regarded as a map into W^ed- 
Let the quadratic Hamiltonian for the infinitesimally hnear map Np^ be Hnn. 
With respect to (pO[), the array for Hnn is 



























^131 

Pe 

^23 

Pe 








^211 

Pe 















































W^)*<(Pe) 























-Id 






Id 














211\* 





^23 

Pe 





where w^^(pe) denotes the symplectic form on W^4d, « = 1,2. Using (|l|), and 



0* 



0^* are certain poly- 



the fact that the projections from the sphtting g*^ — q 

nomials in coad^^ , the bilinear form N^^ block diagonalizes over that splitting 
Denote the blocks by iV^^^° and N^^^. Then, using the variables 



' ' rn 



d X 0u, X 0, 



we have 



1 



1 

2 



(24) 



I now construct the normal form symmetry for iJnii- Let T^^ be the closure 
of exp(R^e), so T^^ is Abelian, is generically a maximal torus of G^^, and T^^ 
naturally acts on both g^^ and g* ^ , with invariant subspaces g^^ and g"^ , i = 
1,2. By (|l3|), Np^ is T^^ -invariant, and so (|2j) is invariant under 2^ restricted 
to the fjP variables alone. After this, it is the last two terms of (|2j) that will 
dictate its further symmetries. Now S-^ also actsQ on g|j__ by 



exp( 



= exp ( 



which is merely a renormalization of the action of T^^ on gj^^ , and by dual- 
ity there is a corresponding action of on g*^. Also there is the following 
symplectic action of S*^ on W^Jd- 



exp(6' )x EE exp ( —dXn^^ {pe) ] x, 9 e 



(25) 



^Define 9^ 



-1 

1 



so that = 50(2) and exp is the usual matrix exponential. 
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By and (jl^), -/Vp^-'^ and Npl^ reverse- intertwine the two actions, so that 



Consequently, the normal form symmetry will be T^^ x S"^ where the first factor 
acts on the variables via the coadjoint action and the second factor acts by 

a ■ (x, jj.^) = {a~^x,afj,), a — exp{9'^). 

The normal form symmetry restricts the quadratic Hamiltonian H^u , as fol- 
lows. There is a basis of the linear space M^J^ = {(2^1 12:2) G K^} so that the 
symplectic form is dx^ A dx^ and the action on W^l^ is either (exp(6''^), x) 1-^ 
exp(6'^)x or (exp(6'^), 1— > exp(— 6''^)a; (i.e. counterclockwise or clockwise ro- 
tation). Using the Ad- invariant metric of g, choose a orthonormal basis for the 
space , so that gj^^ = {(tti, 7r2) € M^} and the action on gj^^ is a rotation in 
the opposite sense as the action on M^ed- The map N^^^ is highly restricted: 
since it reverse-intertwines these two reverse actions, there is a k such that 



1 
1 



or A^2ii ^ ^ 



-1 

1 



(26) 



Equation ( pq ) defines k, so that the map N^^^ has the effect of giving rise to this 
single "coupling constant" and one of the two possible choices in (^6|) . Without 
loss of generality, by the symplectic isomorphisms {xi,X2) '-^ {—xi,—X2) and 



{xi,X2) i—X2,xi) on , K > and the second choice in (26) may be 



discarded, respectively. Then, since N^^ is iS^-symmetric, a constant Ii may be 
obtained from N^^^ such that 

^nii = ^^pf°(/, /) + ^^1(7^1" + ^2^) + + ^2x2). (27) 

According to the above conventions, N'^^^ reverse intertwines reverse-sense 
actions on the {xi,X2) variables and the (7ri,7r2) variables, and the second 
action has the same sense as the action of adj^. Consequently, there are two 
distinct possibilities: the action on the {xi,X2) variables can be counter- 
clockwise or clockwise. To avoid unnecessary signs that would otherwise appear 
later, I will redefine the actions in the following way: if the action on the 
{xi,X2) variables is clockwise then reverse it, and both actions are then coun- 
terclockwise and the action on the (tti, 712) variables is in the same sense as 
the action of ad^^. If on the other hand the action on the (xi, X2) variables 
is counterclockwise then the action of the {7Ti,7T2) variables is clockwise, and is 
to be reversed, so that it has the opposite sense to the action of ad^^. Thus, 
by these conventions, both actions are counterclockwise, N'^^^ is intertwining, 
and the distinctness of the two cases appears as one of two possibilities: the 
action on the (tti, 772) variables might be in the same sense as the action of ad^^ 
(I call this the case, and it occurs when the original action on the {xi,X2) 
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variables is clockwise) or the opposite sense as the action of ad^^ (the "— " case, 
occurring when the original action on the {xi, X2) variables is counterclockwise). 
Although it may appear that same-sense vs. opposite-sense with respect to adj^ 
holds some significance, similar to, for example, the difference between a 1 : 1 
and 1 : —1 resonance as in Kummer [1976, 1978|, this is not so: changing the 
the symplectic form dxi A dx2 to dx2 A dxi will exchange my and "— " 
cases. These two cases are qualitatively identical and the subsequent theory is 
ambidextrous with respect to them. 

Now I will determine the effect on of choosing various subspaces W^cd, 
or equivalently, various operators A commuting with Sp^ in (p^). With respect 
to the decompositions (^8|) and ( [l9| ) let A have the block form 



A = 



^11 ^12 
^21 ^22 



To determine which A commute with Sp^ , note that Sp^ has a block diagonal 
wi m„o - " — ° 



form on W^^^ 



od ^ 0/1 ® 0/^ since (|T8|) and (|19|) are generalized eigenspace 
decompositions for Sp^ ; let the blocks be^^o, S'"Vs'220 , and S^'^^ ,soAS = SA 



becomes 

^21 



^12 
^22 



5110 








^220 








^11 

^21 



^12 
^22 



(28) 



When multiplied, the top left blocks of (H) give A" 5""" = S'^^o^ii^ which is 
possible only if is zero, since S^^'^ and 5"^^*^ have no common spectrum, and 
similarly A^"^ = and A^^ = 0. However, 5^-'^^ and S'^'^^ both have spectrum 
{±Ae}, so one cannot conclude A^'^ = 0, but only that A^^^-in ^ 5'22i^22 gj^^^g 
S^^^ is exactly the hnearization dXH^,^{pe) restricted to W^\^, which by ( p5| ) 
defined the action on W^.^di which is counterclockwise or clockwise, we have 

1 

-1 



Since 5'^^^ is — ad^_^ on q]^^, and ad^^ generates the 5*^ action of opposite sense 



on 

to that generated by S , we have 
By letting 

^22 



1 

-1 



a^i a22 



^225.111 ^ 5-221^22 becomes a^^ = -a^i and a" = 0^2. By (|2^) and the 
duahty M), 





K 

~K 
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so that 





-K 





2ai2 
2ai2 



so by (JlTl), the Hamiltonian generated by the new choice of W^cd has the same 
form as ( p7| ) but with Ii replaced by h + 2na^^. Thus, the entire effect of 
the freedom to choose Wrcd is exactly that the constant Ii may be arbitrarily 
manipulated, as long as k ^ 0. The effective action of the gauge group as it acts 
on the normal form is isomorphic to the additive group of real numbers. 



3.2 The normal form 

The splitting of Tp^P and subsequent analysis of the linearization dXij^^ipe) 
amounts to the construction of a normal form for linearizations of relative equi- 
libria. When equivariantly moved through G^^p^, this normal form can be com- 



bined with th e equivariant isotropic embedding theorem ( Marsden | 1981 |, We- 
instcin [ 1977| ), to provide a G^^ -equivariant normal form for the entire Hamil- 
tonian system near G^^^pe- I now construct this normal form, closely follow- 
ing |Patric| |l995t . 

Let P be the symplectic manifold 



T*G^, X G^e 



and define 



Pe 



e ® Old ® Me 



Let the group G^^ act on P by left translation on the third factor. At the 
zero section of T*G^^, there is the splitting of the tangent space of T*G^^ into 
vertical and horizontal parts, and hence 



'0A.e ©fll ®^fe(GMe). 



(29) 



Together, (poj ) and (P9| ) yield a linear symplectomorphism between Tp^P and 
Tp^P, which extends by G^^-equivariance to a vector bundle isomorphism, say 
A' : TP\{Gfj,^Pe) TP\{G^^Pe), which is symplectic on each fiber. The equiv- 
ariant isotropic embedding theorem then extends A' to an equivariant symplec- 
tomorphism A from a neighborhood of G^^pe to a neighborhood of G^^p. The 
relationship between A' and A is A' = TA on the domain of A'. 

The linearization dXn^^iPe) is in a sense a first order approximation of 
Xh/;^ Sit Pe, and the quadratic Hamiltonian associated to dXH/,^ (pe) is a second 
order approximation to i?^^ at Pe- Together with the map A, this yields a 
second order approximation to the Hamiltonian 7J in a neighborhood of G^^pe- 
Indeed, let be the momentum mapping for the left action of right translations 
(g, h) ^ hg^^, so that 



J«(ag) 
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where is the left invariant vector field generated by Then 

H = Hail — ( J^" — i?rod) + -^rod + ^rsd + ^^rmdj (30) 

where 

1. Hnii : W^l^ X T*Gf,^ R is the extension of (|2|) to r*G^^ by left invari- 
ance, 

2. : W^rcd K is c^^i?^, (Pe)|M^/ed, i = 1, 2, 

3. H,sd ■ ^ M by H^sdigfJ-e) = -(5/ie,Ce>i and 

4. difimd = and d^H^md = on G^^Pe- 

Items (1)~(3) are definitions, while the content is in item (4), which follows since 
the two Hamiltonians 

both have relative equilibria at pe and have been manipulated to have matched 
linearizations (and consequently matched Hessians) at each point of G^^pe- The 
point is that a second order approximation to H may be obtained by dropping 
the remainder term i?rmd of (|30|). 

Now, consider 11 as defined by (pO|). The second term (the one grouped 
within brackets) generates the one parameter group of symplectomorphisms 
generated by <^e through the action (|3l]). Therefore, the first term i/nii and the 
second term Poisson commute. Obviously, then, the first four terms of ( pO[ ) 
pairwise Poisson commute. The effect of the term — J^" is merely to generate 
the flow of the right invariant vector field generated by ^e, and this corresponds 
to the "fast" evolution of the relative equilibrium pe , while the flows of the terms 
-ff^gj and -fffsd cannot cause evolution along G^^pe- So the slow evolution near 
Gfi^Pe is to flrst order dictated by the flow of i/nii- 

As for the symmetries of H, since the Hamiltonian H is G invariant and A' is 
Gyj^-equivariant, the Hamiltonian ( ^o| ) is G^^ -invariant, including the remainder 
term. Also, when Hni\ is extended to be a left invariant function of M^J^j x T*Gij,^ , 
invariance under the normal form action given by ( |l6| ) becomes right invariance, 
i.e., for the part, invariance under the action 

a ■ {x,ag) ^ [a^'^x,{TRa)*ag), (31) 

where R denotes right translation of G^^, and the functions H^^^ and Hisd 
are also invariant under (|3^). However this normal form symmetry does not in 
general extend to the remainder term i?rmd- The action ( |3l| ) is counterclockwise 
or clockwise diagonal rotation on the variables (771,772) and {xi,X2) in a way 
that is a sense-preserving renormalization of the coadjoint action of exp(K^e) 
on (771,772). Since ( ^ ) is the action of the Abelian S^, this action can be 
reversed and thus assumed to be counterclockwise; however, if this is done then 
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the resulting action may acquire the opposite sense to the coadjoint action of 
exp(R5e) on (tti, 7r2). 

In concept the map A is hke a coordinate system near G^^pe] it is only that 
its range is a manifold rather than an open subset of Euclidean space. Through 
its dependence on the splitting of Tp^P, A depends on the choice of W^od, so 
the effect of varying Wred like varying a coordinate system near G^^pe- Thus, 
the freedom to vary W^ed is like a gauge freedom. I will call a particular choice 
of W^ed 3 choice of gauge, and say that, in view of Section (3.1), the gauge 
group of pe is (R, +). As shown in Section (3.1), varying V^-ed itself is a linear 
process, the effect of which on H^n is easily found and described. Yet, through 
the isotropic embedding theorem, A depends on W^od in a nontrivial way. 

This ability to cause very nontrivial "coordinate changes" while easily cal- 
culating the concomitant effect on Hi^w will enter powerfully in the subsequent 
analysis. Particularly, one can imagine choosing a gauge that "simplifies" Hm\, 
making the analysis of iJnii tractable and yielding information about the flow 
of the original Hamiltonian system near G^^pe- On the other hand, a simple 
-ffnii will in general be nongeneric, so that higher order terms in i/rmd might de- 
stroy whatever delicate nongeneric structures arise from -ffnii- In that case, one 
can focus on the identification and perturbation of these nongeneric structures, 
which will occupy smaller regions of phase space as ones attention is restricted 
more nearly to G^^Pf.. 

Summary: The phase space Pdrift on which some Hamiltonian -ffdrift approx- 
imates to first order the evolution of the system, when that system is perturbed 
from pe, is 

Hdriftix, ag) = {g^^ag,ie) + \Nll{ag,ag) -f {g-^ag^Nf^). 

The equations of motion for this systerr^ are 

^idq ^ d-K , dx , 

where 

n = Ce + {N^^tn + Nll^x. 

There is a splitting — {(tti, 1:2) G K^} ® {m*^} such that i?drift becomes, after 
Poisson reduction by left translation ofG^^, 

Hdrm = Ce) + ^iVpf°(/, /) + lh{7Tl^ + ^2^) + ^^^iXl + ^2X2)- (32) 



^ For a Hamiltonian system on the cotangent bundle of a Lie group G of the form ag 1— > 
l/2h{ag, Og) + {ag, X) where ?t is a left invariant metric on G and Jf is a left invariant vector 
field on G, one has the equations of motion 

f2 = /('tt + X(e), TT = g^^ Og , — = Q, — = — coadn tt. 

dt dt 
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Here iV^^" is a constant quadratic form, and Ii and k are constants, and all 
of these are calculable from the nilpotent part of the linearization of the relative 
equilibrium. The Hamiltonian is invariant under counterclockwise diagonal 
rotation on the variables (7ri,7r2) and {xi,X2) and the coadjoint action of T^^ 
on the variables fjP . The gauge group of the approximation is (R, +), and the 
effect of the gauge freedom is to exactly to undetermine Ii . 



4 The mechanical case 



Here I consider the special case G — SO (3) x (5*^)". Standard identifications 
give g = M'^ X R" = 0*- In order that the momentum be nongeneric, assume 
that the 5*0(3) part of the momentum is zero. Reorient the system (i.e. left 
translate it) so that = {(tti, ■7r2, 0)} C so(3). Since the action of adj^ fixes 
flj^^, the SO{3) part of is either parallel or antiparallel to to k G so(3), where 
k = (0,0, 1). I drop the first term of (p2|), since the background motion of the 
relative equilibrium pe itself will not be of interest, after which ( p2[ ) becomes 

-ffdrift = + ^^2^) + ^-^27r3^ + K{niXi + 1^2X2) 



+ TTa ^ aiPi + ^ OijPip-j, 



i=l 



where the variable labeling is 

(7ri,7r2,7r3,pi, . 



and Gi and are constants. 



The whole of (5^)" is ignorable; an Abelian reduction yields the phase space 
and Hamiltonian (which by abuse of notation will have the same names) 



^drift 



X T*SO{3) = {{{xi,x2) eR^Ae SO{3),{ni,TT2,TT3) eS.^)} 



-ffdrift = -^IliT^l^ + 71'2^) 



-h'^i' + ni-KiXi + 7r2a;2) 



(33) 



where a is constant. The normal form symmetry is counterclockwise diagonal 
action of on the variables {xi,X2), (tti, 7r2), and the momentum mapping for 
the normal form symmetry is 



r 



\{XI^+X2^) 



An easy verification (see the footnote on page |l^) gives the equations of motion 
for (H) as 



d 
di 







I7T2TT3 ~ KX2TT3 + 0^2 


71'2 




— ItTiTT^ + KXlTT^ — flTTl 


7I"3 




K{niX2 ~ TT2Xl) 


Xi 




K-K2 


X2 







(34) 
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and 



A- 



dA 
'dt 



/iTTl + KXl 
IlTT2 + KX2 

/27r3 + a 



(35) 



where I = I2 ~ Ii and as usual, for v G 








_^3 






^3 

















To obtain a first impression of ( p3|) one can contrast it with the nonreso- 
nant case. On the one hand, the phase space Pdrift of (|33| ) has dimension 8; 
the symmetry group is 5*0(3) x S^. Hence, the Marsden-Weinstein reduced 
systems are integrable, since they have dimension 8 — 4 — 2 = 2, and energy is 
a conserved quantity. We shall see that homoclinic connections appear in the 
reduced spaces of (33), which are spheres in general, and that singular reduced 
spaces (pinched spheres) also occur for some momentum values. On the other 
hand, the nonresonant case (obtained from (^^ by setting k — and dropping 
the variables xi and X2) has the same symmetry group but dimension 6 phase 
space; it is integrable by its symmetry group alone, without the use of energy. 
Dynamically, the nonresonant case is simpler, since every motion is a relative 
equilibrium and the compact symmetry group rules out homoclinic connections. 



4.1 Full reduction 

The Poisson manifold x R'^ = {(a;,7r)} may be reduced by the normal 
form symmetry by use of the Hopf variables (see, for example, Cushman ancj 



Bates [1997 1, page 14; also page 407 has a summary of references about singular 



reduction): 

Wl = 2(xi7r2 - X27ri), W2 = 2{xi'Kl + X2TT2), 

W3 = W + 2^2^) - (vri^ + vra^), W4 = (xx^ + X2^) + ij^i + t^2^). 

This map is a quotient map for the normal form symmetry and it has the 
semialgebraic image given by the subset of — {(wi, W2, W3, W4, 7r3)} satisfying 

•W\' + W2^ + Wz' — W4^ =0, W4 > 0. 

The symplectic reduced spaces are the level sets of the two Casimirs j\ and j2 
given by 

jl EE (7ri2 + 7r22+7r32)^ ^ ^7r32 + ^(w4-W3)^ , (36) 

h = ^"^ = -^i^l^ + ^2^) - 713 = -^{W3 + Wi) - 713. (37) 
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Given ji and j2, ( |36[ ) and ( pT\ ) can be used to eliminate W3 and W4, thereby 
obtaining the symplectic reduced spaces as the subsets of = {{w 1,-0)2,^3)} 
given by 

w^^ +W2^ = 8(713^ -ji^)(Tr-i+j2) (38) 
W4 - (il' - 2j2) - (vrg^ + 2^3) > 0. (39) 

Various symplectic reduced spaces are obtained by fixing various values of 
ji and j2 , but the values of ji and j2 are not arbitrary, since ji > and 

il = (tTi' + TTa' + TT3^)i > -TT3 = j2 + ^{xi^ + X2^) > 32- 

Moreover, the singularities in the symplectic reduced spaces occur where xi = 
X2 — TTi ^ 712 = 0, and here ii = Ivral = ±i2. Equation ( ^ ) has solutions only 
over the intervals where the cubic on the right hand side is nonnegative, and 
moreover, by (|37|), tt^ > —j2- Putting all this together gives the bifurcation 
diagram shown in Figure (|l|). As shown there, the symplectic reduced spaces 

are the surfaces of revolution of (8(7r3^ — ii')(7r3 +i2)) ^ over the finite interval 
where the cubic is positive. Thus, these spaces are points, spheres, or, in the 
case that i2 = — ii, they are topological spheres with one conical singularity. 

In passing, I note that the symplectic volumes of the — i2 = ii > reduced 
phase spaces are infinite, and of course the infinity is concentrated at the singu- 
larity. In fact, by regarding two functions f{wi, W2) and .g(wi, ^2) a-s functions 
on one of those reduced phase spaces, lifting those functions to the Poisson 
phase space {(a;,7r)} = x so(3), and then calculating the Poisson bracket, 
one verifies that 



{f,9} = -4(713' + 2i27r3 +ii 



/ df dg df dg 



\dw2 dwi dwi dw2 



Consequently, the symplectic form on the reduced spaces in the coordinates 
(wi, W2) is 

4(7r3^ + 2j27r3 + Ji ) 

Switching to polar coordinates wi — rcosO and W2 = rsinfl, the symplectic 
volume of a small circle of radius S about the origin is 

T TT 

Aa = 27r / ^ dr ^ - I '■ -rr dr 



/o 4(7r32 + 2i27r3+ii') 2 J„ (7:3-31) 

and upon substituting 

= 8(713' -ii')(7r3 -ji) = 8(7r3 +ii)(7r3 -ii) 
that symplectic volume becomes 

A, ^47: r^^^±^dr. 
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Figure 1: The bifurcation diagram for the symplectic reduced spaces. Clockwise from the 
top, the reduced spaces are points, spheres, pinched spheres, again spheres and finally points. 
The finite intervals uihe re t he cubic is positive is the thicker black line on the Tra axis. One 
verifies that Inequality (Isdj is respected by showing that the concave down quadratic W4 = 
— ("■3^ + 2TT'i) — ji'^ + 2j2 is positive on those intervals by checking the it is positive at the 
various endpoints of the intervals. 
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Since the numerator tends to 2ji as r becomes small, Ag is infinite. 

The Hamiltonian (^) on the reduced spaces is, with the aid of the Hopf 
variables and (|3^), 

1 Ik 

-ffdrift = jhiWi - W3) + -hTTs^ + -^W2 + a-Ks 

= ^1^3^ + + a.TT3 + ^hji^ 

Obviously this expression is simplified if Ii ~ I2, since then I = I2 — Ii =0. 
This choice of /i gives the first part of ( |33| ) the same form as the kinetic energy 
of a spherical ball moving in 3-space, so I call this choice the spherical gauge. 

In the spherical gauge the flow lines of the reduced systems are trivial to 
determine: they are the intersection of the planes {k/2)w2 + otts — constant 
with the surfaces of revolution in Figure ( [4.1| ). In the case of a nonsingular 
reduced space of nonzero dimension the flow is that is that of the flow on a 
2-sphere with two stable equilibrium points. For the singular reduced spaces 
— j2 = ji > the type of flow depends on whether or not the plane cuts through 
the cone at the tts = ji singularity of u;2^ = 8(713 + ji)(7r3 — ji)^. That cone is 

1^2! = 4^71(^1 - TTg) ^ TTg = ji - --^|u;2|, 

4VJi 

and so there is an unstable equilibrium at the singularity if the line {k/2)w2 + 
aiT^ — aji passes through that cone, which is when < 4k^ji, and a stable 
equilibrium otherwise. The former case is the flow on a 2-spherc with two stable 
equilibrium points and a single homoclinic connection to a singular point, while 
in the latter case this flow has two stable equilibrium points, one of which resides 
at a singularity. 

For the remainder of this section I will assume the spherical gauge. 

The equilibria of (|^) correspond to relative equilibria for just the action 
of SO{3), since the space {{tt,x)} is the Poisson phase space of Pdrift reduced 
by the SO{3) action. These have a greater significance than relative equilibria 
which also use the normal form symmetry, since the full unapproximated 
Hamiltonian has the 5*0(3) symmetry but it is only the truncated approximation 
which has the 5*0(3) symmetry and the normal form symmetry. As well, it 
is expedient to separate the equilibria of ( ^ ) at the outset of the analysis. 



Setting the right side of (34) to zero immediately gives tti = 7r2 = and then 



XiTTs = X2Tf3 = 0, so there are the following solutions: 



X = 0, TT (40) 

a; = 0, TTi TTa = 0, 713^0 (41) 
a; ^ 0, TT = (42) 



The set of relative equilibria given by (|4^) contains relative equilibria equiva- 
lent under the 5^ normal form symmetry. The 50(3) generator corresponding 
to (E2h is, by substitution into (pd), {nxi, kx2,0). The zero section of T*50(3) 
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(i.e. the set tt = 0) corresponds to perturbations of the original relative equi- 
librium such that the perturbation has zero total angular momentum, and 
in the nonresonant case such perturbations would imply no drift and no drift 
has been numerically observed for such. For the resonant case, however, one 
expects regular rotation around the axis (kxi , kx2 , 0) , which is perpendicular to 
the generator ^e- Nonzero x corresponds to some "reduced excitation". Conse- 
quently, at zero total angular momentum the system drifts so that the generator 
of the relative eguilibrium moves along a fixed great circle at a rate dictated by 
a reduced excitation. 

The relative equilibria ( |4^ ) are the relative equilibria with nongcncric mo- 
menta for the system (|3^). Consequently, the methods of Patrick | 1995 | might 
be applied to them: the drifting relative equilibria could themselves drift. I will 
not, however, pursue this aspect here. 

To obtain a list of nonequivalent relative equilibria (i.e. a list of relative 
equilibria, no two of which arc in the same SO{3) x orbit), assume yl = Id 
and discard Equation (|35|), and, using the normal form symmetry, assume 
TT2 = and TTi > 0, and when tti = assume that X2 = 0. Equating the first 
of (M) with the infinitesimal generator of Se € M gives 



(43) 



Putting = into the fourth component gives Si,X2 ~ 0. Now = corre- 
sponds to the 5*0 (3) relative equilibria (|^-(||), and in ^) and (|l|), we have 
X2 — 0, while in ( [42| ) we can assume X2 = and xi > by the normal form 
symmetry. So if Sg = then we assume X2 = 0, while if Se ^ then X2 = 
anyway by SeX2 = 0. Putting X2 — tt2 = into ( ^3|) then gives 



-KX27r3 + a-TT2 






KXiH^ — flTTi 






k{tIiX2 - 1^2X1) 







KTT2 






— KTTl 




SeXi 



and these equations are easily solved to obtain the following list of nonequivalent 
relative equilibria: 



TTi = 7r2 = xi = X2 = 0, r]e = (a + .Se)k, TTs e M; (44) 

TTl(aXl — KTTl) KTTl 

KXl^ ' Xl ' (45) 

rje = KXii + {a + Se)k, tti € M, > 0. 



7r2 = X2 = 0, TTs 



Here the 50(3) generators r]e of these relative equilibria have been determined 
by comparison of ( |35| ) and the infinitesimal generator of the 5*0(3) action of 
-fdrift at the relative equilibrium. In (|4|) rj^ and Se are not unique due to the 
presence of isotropy in the phase space Pdrift- For later use, the characteristic 
polynomials of the linearizations of these relative equilibria are, respectively, 

Pi = x^{x^ + \rie\'^){x^ + (a^ + 2ase + 2s^ - 2t:3K^)x^ + {sj + as^ + ^3^^)^) 
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and 

P2 = -^x'^ix'^ + \7-ie\'^)(xi'^x'^ + K^Xl'^ + {axi - 2K7ri)^). 
Xi'^ 

For an understanding of the relative equilibria and the flow in general it 
is important to understand which relative equilibria reside on which reduced 
spaces. For (|4|), ji = IttsI and j2 — tt^, so for tts < one has ji ~ j2, which 
corresponds to pointlike reduced spaces, while for tts > one has ji = —j2, 
which corresponds to the singular reduced spaces, and the relative equilibria 
occupy the singularities. As for (|45|), substituting ( |4^ ) into ji and j2 and elim- 
inating TTi gives 

(3xi4+8j2Xl2+4j2'-4ji2)2K2 

+ Aa^xi^ix,^ + 2j2 - 2j,)ix,^ + 2j2 + 2ji) = (46) 

so it is a matter of solving this quartic in xi^ for its positive roots. The pointlike 
reduced spaces corresponding to ji — are occupied by the relative equilibria 
of (45) after putting tti = 0, corresponding to the solutions ji = 0, a;i = ^—2j2 
of @. 



4.2 Partial reduction by the normal form symmetry 

In this section I symplectically reduce the Hamiltonian system by the normal 
form symmetry at a value, say a of the momentum mapping J"^. My objective 
is an interpretation of the reduced system via the equations of motion on the 
reduced space. 

Use the notation = xi'^ + xi^. The cr-level of the momentum map J'^^ is 

1 2 

TT-? = —a r . 

2 

Below, when x appears in the context of a vector in M"^, it is as (xi, X2,0). The 
map 

^'(A, TT, x) = {y, a, z) = (^k, A(7r - (tt • k)k), Ax) (47) 

is a projection of the phase space T*SO{i) x to the Whitney direct sum 
j'*g2 TS"^. The restriction of this map to the a level set of J"^ is clearly a 
quotient map for the normal form action, and thus the symplectic reduced space 
is r*S'2 © TS"^. The action of £'0(3) on T*S^ ® TS"^ becomes 

A - {y,a,z) = {Ay,Aa,Az) 

To find the reduced vector field at (y, a, z) e T*S'^ ® TS"^, first evaluate the 
original Hamiltonian vector field, namely 

f^EE/i7r + /((7r-k) + a)k, A^^^ = f^^, 

at 
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= TT = TT X + ak, 



da: 



= X = — /ck X TT 



(48) 



at the point 



■K = A ^a- (T+-r^ k. x = A^z 



(49) 



where ^ is chosen so that Ak — y, and then apply the derivative of (47), which 
is 



T-^{A,Tr,x,Q,n,x) 

= -^{Aexpin^t),?: + tTr,x + tx) 

= A ■ (n X k,n X (tt - (tt ■ k)k) + tt - (tt • k)k, x + n x x) 



(50) 



Thus it is a matter of substituting ( p8| ) and ( p^ ) into (^ . Without care the cal- 
culation can be onerous; however, it is straightforward, and yields the equations 
of motion 



— = [ha + Kz) X y, 
at 

da 

— — —TT^llia + KZ) X y + KZ X a, 

dz 

— = (na — (/27r3 + a)z) x y + Iia x z, 
dt 



(51) 
(52) 
(53) 



where tt^ stands for —{a + |zp/2). To these equations must be added the 
constraints \y\ = \, a ■ y — Q and z ■ y — Q. 

I want to impose the viewpoint that there is a "particle" at y having direc- 
tion z. To do this, I use to replace a with dy/dt in ( |5^ ) and ( p3| ) while 
writing equation (|5^) as a second order equation in y. I also use the standard 
Levi-Cevita connection of S'^ for the time derivatives. Thus, Equation ( p2| ) 
becomes 



df^ 



dy 



{Iia + Kz) X — + [ h— — h K 



dt 



da dz 



dt dt 



X y 



= (/la + Kz) X ((/la + Kz) X y^ 

+ /i7r3(/ia -I- kz) - K^Ka — {hTTs + a)z) 



dy 
dt 



dy 

y + iiT^sy X -77 - 

dt 



y X ^ - kz] ~ (72773 -I- a)z 



and similarly with Equation (pS^, so that 



W^y /i V3 - k2 dy 



dt^ 



^2773 + a] z 



Vz K dy f K^ . 



(54) 
(55) 
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Some aspects of the particle become apparent by replacing a with dy/dt in 
the total energy and angular momentum. Using the total energy is 



2/1 



-7r3 + flTTa + Kz ■ a 
\Iia + Kz\ + -TTs^ + a-TTa - 



2/1' 



(56) 



1 


dy 


' 1 / 


2h 


dt 





7:1- -;-+ I2T^3 



h 



SO the internal energy of the particle is 



1 

2^2 



The total angular momentum, or the 50(3)-momentum-mapping is, using (p^, 
J = Att = a — TTsk = — y x — —z — vrsy, (57) 



1 dy K 
TTsk^ —y X — - —z 
li dt li 



and below I replace z with the tangent to the sphere part of the angular mo- 
mentum, which by (^^ is 

r - 

L = -—z, 

while the total angular momentum attributed to the particle must be L + Tr^y. 



By the first terms of both ( |56| ) and (57), the particle should be viewed as having 
mass 

From Jackson |1975] the classical nonrelativistic equations of motion on \y\ = 
1 of a particle of mass m, charge Q, and gyromagnetic ratio F, and magnetic 
moment m = FL, are 



m- 



dt^ 



— B X V + V(m • B) - (V • B)m, 

c 



= I + TLx B, 



(58) 
(59) 



where c is the velocity of light and v is the particle's velocity. Here 2 is compen- 
sates for the inductive forces that would, were m to be constant, be required to 
maintain the constant current that generates m itself. One can fit equations ( p4| ) 
and (55) to these equations as follows. If one takes B — y, comparison of the 
first terms of (pq) and (m) gives the inductive term 



(60) 

if the gyromag- 
(61) 



_if_dy 
dt ' 

while the second term of (^5|) matches the second term of 
netic ratio is 

2 



K 
77 
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Matching the second term of (|5J) with the third term of (|5^) imphes F is «;^/3Ji 
times the left side of (|6l|). So choose a gauge so that Ii — Then matching 



the first term of (54) with the first term of ( pq ) yields 

Q 

C 

Finally, I note that the inductive term ( |60| ) is exactly what is required to balance 
the rate of change of internal energy with the work done on the particle through 
its magnetic moment via the last term of (|5^). 



5 Example: Two coupled rods 



In this section I illustrate and verify the above theory using numerical sim- 
ulations of the system of two axially symmetric rods which are joined by a 
frictionless ball-and-socket joint. In the center of mass frame this system can 
be cast as a geodesic flow with configuration space S'0(3)^ = {(^11,^2)} and 
kinetic energy metric 



^2' 



J{A) 



17i 

^2 



1 


















1 


















a 


















1 


















1 




















a 



where A = [A^^] = A1A2, ili and are body referenced (left translation) 
angular velocities, and 



J{A) = 



Here 0<Q!<2 — 2/3isa parameter which increases with the diameter of the 
rods and < /? < 1 is a parameter which measures the degree of coupling 
between the rods; /3 is zero if the joint lies at the mutual centers of mass of the 
rods. The continuous symmetry group of this system is 5*0(3) x (S"^)^, which 
acts on the configuration space 50(3)^ by 

(B, 01, 02) • (Ai, ^2) = {BAi exp(-0ik'^), 5^2 exp(-02k^)). 

All relative equilibria for this system are explicitly known, as are their formal 
stability; a complete list can be found in Patrick| | |1991 1 . The relative equilibria 
of interest here are the phase space points (Id, exp(0j'^), Vli, ^,2) parameterized 
by ti, t2, e R such that ti 7^ 0, i2 7^ 0, < < tt, where 



ill = tii 



k, ^2 



t2i 
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and where 



ati sin 6 

The corresponding generators (il, cti, (72) E so(3) x are given by 

CTl - , (72 - , - Ell H ^T^j^ K. 

By calculating the linearizations, one sees that a 1:1 group reduced resonance 
with zero total-angular-momentum can be arranged by setting 

2t t 

P = 3 ^ % , 3ti^ + ti%^{Acos^e - 10) + 3^2^ = 0, 

tl" + t2 

and for the work presented in this section I have chosen to perturb the particular 
relative equilibrium obtained by setting 

TT 1 + 75 1 2 

0=-,t,^l,t2^^-,a=-,(3^^. (62) 

5.1 Calculations for the coupled rod system 

To compare the dynamics of the drift system to the dynamics of the coupled 
rod system near the relative equilibrium (^), it is necessary to calculate the 
splitting or (H) in this special case. Here is one general way, inspired by the 
proof of Theorem 3.1.19 of Abraham and MarsdcE 1978| : 



(a) Calculate the semisimple part Sp^ and the nilpotent part Np^ of the lin- 
earization cIXh^^ (Pe)- 

(b) Since Sp^ is semisimple it has a basis of eigenvectors, some complex and 
some real, and the complex eigenvectors may be grouped into complex 
conjugate pairs. Taking the real and imaginary parts of one eigenvector 
in each pair, and then including the real eigenvectors yields a basis B of 
Tp^P in which Sp^ is skew symmetric. 

(c) Let the matrix of uj{pe) with respect to the basis B be —W. Since S is 
infinitesimally symplectic and skew, = S*W + WS = —SW + WS, so 
S and W commute. By going to a basis of eigenvectors of the positive 
matrix — W'^, find a symmetric, positive square-root B of —W'^ such that 
B and S commute, and set J = WB~^. Using the basis B, regard J as 
its corresponding operator on Tp^P, and regard B as the bilinear form on 
Tp^P corresponding to {x,y) ^ x*By. Then 

uj{pe){v,w) — B{Jv,w), = — Id, 

and J commutes with Sp^ . 
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(d) Set Z = J{g^J, so Z is a Sp^ invariant Lagrangian complement to g^^ in 
g^^ © Z. Choose an Ad-invariant complement b to g^^, and set V^-ed = 
(0M=-Pe (B Z (B bpe)'^^- Then W^ed is symplectic, Sp^ invariant and is con- 
tained in kerdJ{pe), since Vl^ed C gp"^. The subspace W^cd complements 
g^^ in kerdJ(pe) since g^^ C W^'^^ and W^cd is symplectic. 

In the case of two coupled rods, further refinement of the splitting of W^cd to 
resonant and nonresonant parts is not necessary, since the reduced spaces for 
zero total-angular-momentum are two dimensional; there is no "nonresonant" 
part of the phase space. Final adjustments to achieve the spherical gauge are 
easily arranged using ( [T^ ) and ( [l7| ) . The end result of all this is a basis of the 
tangent space at the relative equilibrium which reflects the splitting (H) , a basis 
of the Lie algebra so(3) x R'^, reflecting the splitting ([l^) into resonant and 
nonresonant parts, the information 

K = .9115064, h^h^ 4.321619, (63) 

and the information that (^2|) is a "+" type relative equilibrium. 

Further to calculating the splitting or (H) , means must be found to trans- 
late initial conditions of the drift system to perturbations (i.e. initial conditions 
near (^)) of the coupled rod system. For this a map was used, say ^P, of the tan- 
gent space Tp^P into phase space P such that Tp^'^ is the identity at p^- The 
map \E' only injected initial conditions of the drift system into perturbations 
of (^2|) to first order; actual correspondence could not be easily achieved due 
to the implicit nature of the coordinates provided by the isotropic embedding 
theorem. Moreover, did not match the momentum of the drift system exactly 
to a momentum perturbation of the actual system, while practice indicated that 
an exact match of these momenta is important. For the comparisons just be- 
low, given an initial condition of the drift system with a particular momentum, 
an initial condition of the coupled rod system with matching momentum was 
obtained as ^' of another nearby initial condition of the drift system. That 
nearby initial condition of the drift system was obtained by a slight (second 
order) iterative refinement of the original initial condition of the drift system. 

By the conventions of Section 4, and since the relative equilibrium (|6^ ) is 
of type, the vector (tti, 7r2, tts) = k is parallel to the normalized rotation 
vector of the coupled rod system. Also, the quantity A e 5*0(3) of the drift 
system corresponds to the drifted orientation of the coupled rod system. Thus, 
the quantity A\s. of the drift system corresponds to the drifting rotation vector of 
the coupled rod system. Now, the rotation vector of the relative equilibrium ( |62| ) 
is a constant linear sum of the locations of the rods, exactly because the rods 
when in the relative equilibrium rotate around that vector: if one defines 

^ (tiAik-ta^ak) 



sin 6 



then T = at the relative equilibrium (|6^) and also at any reorientation of ( |62|) . 
Thus, for a perturbation of (p2), r approximates the drifting rotation vector and 
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Perturbation 



0.04 



Figure 2: Perturbations within zero-total-angular momentum. Left: the paths traced out by f, 
resulting from 6 initial conditions within zero total-angular-momentum. Right: the rotation 
rates of t vs. the magnitude of the perturbation \x\. 



the prediction of the preceding theory is that 



The numerical verifications I have undertaken consist of predictions of the drift 
system for the motion of Ak and the comparison of these predictions with 
evolution of the the unit vector f calculated via a symplectic integration of the 
coupled rod system. The particular algorithm used was a (implicit) Riemannian 
leapfrog algorithm (Leimkuhler and Patrick 1996|). 



5.2 Numerical results 

5.2.1 First comparison: zero total-angular-momentum 

By the results in Section (4.1), the part of the drift phase space corresponding to 
perturbations with zero total-angular-momentum is occupied by 5*0(3) relative 
equilibria, and the motion of Ak from initial condition tt = and x arbitrary is 
that of uniform rotation at angular frequency k\x\ along a great circle through 
k and perpendicular to (a;i,a;2,0). Translated to the coupled rod system, the 
prediction is that f undergoes regular rotation in a great circle through f(0), the 
particular great circle regularly rotating as x is rotated. The left of Figure (||) 
shows the motion of r as a; is varied from to through 57r/6 radians by 
increments of 7r/6, so the coupled rod system conforms to this prediction. As 
the magnitude of x is varied the great circles should be traced out at an angular 
rotation rate k\x\, so k can be determined by plotting that rotation rate against 
|a;|; this is done on the right of Figure (||). The rotation rates fit well to the 
curve 



Rate = .9284|a;| + .7750|xp, 
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so from the simulation k is .9284, which is 1.8% off the calculated value of 
.9115064 already displayed in (|63|). 



5.2.2 Second comparison: a stable relative equilibrium of the drift 

system 

For another comparison, I examined the coupled rod system for perturbations 
of (|62h corresponding to being near one of the stable relative equilibria in the 
list (|45|). For perturbations of ( |62|) corresponding to the values 

TTi = .001, xi = .04 (64) 

in ( ^5| ) the predictions of the drift system are that r moves as small periodic 
oscillation of the full reduction of the drift system superimposed on the periodic 
motion exp(77^t)k, where rje is the 50(3) generator in the list (^5|). The two 
predicted frequencies are the linearized frequency of the corresponding equilib- 
rium on the reduced space of the drift system and the rotation frequency of rj^ , 
which are by substitution of (|6j) into (^ and/or (^6|), respectively, 

.04477, .05837. (65) 

Again by substitution of (Q) into (p5[), the motion of r should be such that its 
projection onto the unit vector ?7e/|77e| is near —.5300. 

Figure (||) shows the results of a simulation of the coupled rod system for 



the perturbation of (62) corresponding to the values (|64D. In the simulation the 
normalized rotation vector f moved on the unit sphere as shown in the top left 
of the Figure; the sense of the rotation is clockwise as seen from the top, so 
that rjf, is pointing down and away from you as you look upon the Figure. The 
height of the rotation vector has been graphed in the top right of Figure (^, 
and is visibly quasiperiodic with about two frequencies. The measured average 
height of the motion was .4970, compared with the prediction of .5300. In the 
middle of Figure (|^) is the power spectrum of one of the horizontal components 
of T. As is shown in the table immediately below the power spectrum, the power 
spectrum is consistent with the 3 fundamental angular frequencies 

.04449, .05791, .1687, (66) 

all other peaks being harmonics of these. The first two of these are by far 
the largest peaks (note the vertical logarithmic scale in the power spectrum) 
and they agree with the the predicted frequencies (|65|). The third observed 
frequency in (|6^) corresponds to a far less prominent peak and is likely a fre- 
quency associated to the higher order terms that have been truncated in the 
drift approximation. 

To give an idea of the speed of the drift, on the short simulation correspond- 
ing to the top left and right of Figure (|^), the system rotated about 300 times 
while f rotated about 7 times. Thus the drift was about 50 times slower than 
the original relative equilibrium. Due to this, these kinds of simulations can 
be long time: over the duration of the simulation that generated the power 
spectrum in Figure (0) the rods rotated about 10,000 times. 



32 







-4- 







A B 






c 














.5 


Power 


Frcq. 


Frcq. 


z\ 






Power 


Frcq. 


Frcq. 


zi 




23 


1.195 


.04449 


.04449 


1 








8.766 


.2266 


.2266 





1 


1 


2.935 


.05791 


.05791 





1 





8.848 


.1822 


.1822 


-1 


1 


1 


3.485 


.01342 


.01342 


-1 


1 





8.874 


.1549 


.1553 


1 


-1 


1 


4.871 


.1024 


.1024 


1 


1 





8.893 


.1607 


.1603 


1 


2 





5.812 


.1687 


.1687 








1 


10.3 


.1296 


.1292 


-1 


3 





6.146 


.1239 


.1243 


-1 





1 


10.58 


.1741 


.1737 





3 





6.485 


.07133 


.07133 


-1 


2 





10.63 


.2711 


.2711 


1 


1 


1 


6.817 


.1162 


.1158 





2 





12.25 


.2401 


.2401 


-1 


2 


1 


7.102 


.2132 


.2132 


1 





1 


12.36 


.2846 


.2846 





2 


1 


8.24 


.06596 


.06634 


-1 


-1 


1 


12.91 


.3371 


.3375 








2 


8.742 


.1104 


.1108 





-1 


1 


12.98 


.2926 


.293 


-1 





2 



Figure 3: Some results of the simulation of the coupled rod system for a perturbation corre- 
sponding to a elliptic relative equilibrium of the drift system. Top left: the evolution of the 
rotation vector projected onto the unit sphere. Top right: the time evolution of the height 
corresponding to the top left. Middle: the power spectrum of a signal derived from the drifting 
motion; note the vertical logarithmic scale. Bottom table: the peaks in the power spectrum, 
sorted by power. The second column tabulates the frequencies of the peaks and the third col- 
umn tabulates the harmonics z\A + Z2B + Z3C of the frequencies A, B and C indicated on 
the frequency axis of the power spectrum. 
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5.2.3 Third comparison: a singularity-induced phase jump. 

As shown in Section (4.2), the singularities of the ~ji — j2 reduced phase 
spaces of the drift system are occupied by an unstable equilibria corresponding 
to S0{3) relative equilibria of the drift system. A single homoclinic orbit, trac- 
ing out the intersection of each reduced space with the plane W2 = 0, emanates 
from every such equilibrium. These homoclinic orbits are interesting features 
of the drift system, in part because they provide an opportunity to investigate 
dynamics near singularities of reduced spaces. I will show here that, due to the 
singularity of the reduced space, the drift system suffers a jump in its reconstruc- 
tion phase as initial conditions traverse the homoclinic orbits, and (numerically) 
that this feature of the drift system persists to the coupled rod system. The re- 
construction phase jump shows the presence of Hamiltonian monodromy in the 
completely integrable drift system; for more information on Hamiltonian mon- 



the same reference, as well as Bates [1991|. 



odromy see Cushman and Bates |1997], page 175, the summary on page 403 of 



So choose ji small and consider perturbations of the relative equilibrium 
TT = jik, X = 0, A — Id, 



where by abuse of notation ji serves as the constant value of the Casimir (g6|) 
of the same name. This relative equilibrium corresponds to the equilibrium 

Wl = W2 = 0, TT-i = ji (67) 

on the 5*0(3) x reduced space. Near (|^ the motion on the reduced space 
of the drift system is periodic along the intersection of the W2 = h plane where 
h is small: for a long time the system remains near the singularity, and then 
moves off, passes near tts = ji and then returns to the singularity, similar to, 
for example, the motion of an inverted pendulum. One reconstruction gives 
the motion on the phase space {(a;, tt)} and, since the motion is periodic on the 
5*0(3) X reduced space, there is a well defined reconstruction phase. On 
the phase space {(tt, x)}, the only points on the level sets of ji and j2 that map 
to ( |67[ ) are 

TT = jik, x = 0, (68) 

since tt^ = ji = —j2 implies by (|36| ) and (|37| ) that vri = tt2 = xi = X2 =0. 
Thus, when the 50(3) x reduced system is near (|67|), the reduced system with 



phase space {(a:,7r)} is near (68). An 50(3) reconstruction gives the motion 
on the phase space M x T*SO{3) = {{x,A,Tr)}, starting, say, at ^4 = Id. By 
conservation of the momentum J = Att, the variable A is nearly a rotation about 
k whenever the reduced system is near (|68|). Thus, the motion of the point Ak. 
on the unit sphere is this: for a long time Ak remains near k, then moves off, 
passes near — k, and then returns near to k. This motion repeats but rotated 
( with respect to the previous excursion from k) by a excursion-independent angle 
about k. This rotation is what I mean by the SO{3) reconstruction phase shift of 
the motion. 
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Figure 4: Left: the motion off for a perturbation corresponding to being near the homoclinic 
orbit of the reduced space of the drift system. The vector f begins near the top of the sphere and 
moves off to the left, comes around the sphere, makes an approximate Txjl turn, then moves 
away from you, and goes around the sphere again. On the right: the SO{3) reconstruction 
phase of the coupled rod system as initial conditions are rotated around the singularity of the 
reduced phase space of the drift system. 



By simulation of the coupled rod system I have verified that the gross details 
of the motion of Ak in the drift system also occur in the coupled rod system 
(see the left of Figure (^). This situation is robust; no particular care is re- 
quired, in choosing the coupled rod system's initial conditions to evoke these 
kinds of motions. The robust nature under addition of higher order terms of 
the homoclinic orbit is expected through averaging theorems such as the one 



in Guchcnheimer and Holmes |1983|, page 168. 

It is the behavior of the reconstruction phase of the drift system near to 
TT — that is relevant for predictions of the motion of the coupled rod system 
near to its resonant relative equilibrium. One way to examine this behavior is to 
numerically integrate initial conditions of the drift system that are reconstruc- 
tions of initial conditions starting at a small circle surrounding the singularity of 
the iS'0(3) X reduced space. For the reconstruction phase shift the result 
is a constant phase shift of ±7r/2 with a jump of tt as the homoclinic orbit is 
traversed. The corresponding phase shift of in fact becomes undefined as 
the perturbation vanishes. The reason for this is that the equation of motion 
for A, namely Equation (|35|), becomes, near to (|68|), the equation 

while A{t)'k is not exactly k due the the presence of the perturbation itself. 
Consequently the ^k picks up a rotation of angular frequency /2 ji + a during 
its long visit of (|6^). As the perturbation vanishes this long visit becomes an 
eternity and this undefines the 5*0(3) reconstruction phase shift as an asymp- 
totic effect. However, viewed from a frame that corotates with the same angular 
frequency, namely l2ji + a, one can expect a well defined phase shift. I have 
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numerically verified this: from the corotating frame the SO{3) reconstruction 
phase shift of the drift system is nearly ±7r/2 with a jump of tt as the homo- 
clinic orbit is traversed. This compares favorably with the SO (3) reconstruction 
phase shift of the coupled rod system shown in right of Figure (||). 

Finally, so that the comparison of these phases of the drift system and the 
coupled rod system are not entirely numerical, I give a calculation showing that, 
asymptotically as the homoclinic orbit is approached, these phases arise mostly 
from the singularity of the reduced space. The reconstruction phase of the 
drift system may be calculated by a slight modification (to allow reparame- 



terization) of the usual reconstruction method found in Abraham and Mars 
den il97S|] : Generally, suppose H is a Hamiltonian on a symplectic phase space 



P and that a curve c^{s) is a reparameterization an evolution c^(t) on some 
Marsden-Weinstien reduced phase space — J^^{^)/G^. Choose a curve 
d{s) S J~^(/i) such that c — ir^ o d, where tt^ : J^^{p) — > Pp, is the quotient 
projection. Then there is a unique smooth curve ^(s) S and a unique smooth 
function a{s) such that 

XH{d{s)) = f(sMs) + a{s)d'is). (69) 

By differentiation, the curve g(s(t))(i(s(i)) satisfies Hamilton's equations on P 
if 

-idg 1 <-/ \ ^ 



ds a{s) ds a{s) 

Perfectly obvious generalizations hold if d is just a 1-manifold in (/i) covering 
the image of c^. Particularly, by setting tt2 = 0, d can be the subset of Pdrift 
defined by 

7^1 + ^2 = jh =0, Xl^ X2 = , (70) 

wl + h^ = 8i7Tl-f,){n,-j,), (71) 
in which case one calculates using (|4|), (|69|), ( |70| ) and ( [7l| ) that 
d9 —hn^ dTTs kwi 



dT^3 (ji^ - Tra^)^! ' dt 2 



(72) 



where 9 is the reconstruction phase as the counterclockwise angle of the 
vector (TTi{t),TT2{t)) from the tt2 — axis. 

Let ri{h) < r2{h) < r3(ft) be the roots of the cubic 8(713 —ji)^(7r3 + ji) — h"^. 
By standard perturbation arguments 



1 



r2^n--r^h + 0{h'), 
4vJi 
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r3 = ji + 



h + 0{h^). 



8(7^1 - jf)i7T3 - ji) is 



The 5*^ reconstruction phase over the curve wl + h? 
calculated as follows. Let the reduced system start at wi — 0, tt^ — r2 at time 
t = to-, move to wi = 0, tts = ri at time t = ti and then complete its periodic 
orbit by moving back to wi = 0, tts = r2 at time ^2- Then the contribution 
of ([t^ to phase shift over the interval [tc^i] is 



0(ti)-0(io) = 



dt 



dt = 



7r3(to) 



(i0 dt 
dt diT-i 



■ dns 



where the positive square root must be used when solving for wi in (|71), since 
TTs must immediately decrease after time t = to, and by the second of (72|), wi 
is positive over the interval [<o,ii]- The contribution of (|7|) over the interval 
[ti,t2] is identical (use the negative root here). Thus, the total phase shift 
cj} = (t){t2) — 4'{to) over the loop as the loop approaches the homoclinic orbit 
through /i > is 



-2 lim 



(ji^ - 7i"3^)v/8(i'3 - ri){r2 - 7r3)(r3 - tts) 



: dlT^i 



(73) 



For small h the integrand of ( |73| ) is small away from its two singularities at 
TTs — Ti and TTa = r2 . I begin with the left singularity at tts = ri , so I calculate 



hi = lim 



1 (ji^ - 7i'32)^(7r3 - ri)(r2 - 7r3)(r3 - tts^ 



Elementary estimates show that zero error is made as /i — > by the replacement 
of 

71'3 



with its limit as h 
lim 



(ii - 7i"3)^(7r3 - r2)(7r3 - rs) 
of its evaluation at "n^ — ri, which is 
ri -ji 



h^o+ (ji _ ri),/(r2 -ri)(r3-ri) 2ji{y/Tny 



-1 
4ji' 



Thus 



hi = lim 



-1 



/i^0+ 4ji\/8 7ri (ji + TTs) V's - ri 



dn-i 



= lim 1= ,. 

h^a+ 4jiV8 VJi+ri 

h 2 vr 

~4ji\/8 ,n^2 

TT 



arctan 



7r3=0 
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The right singularity of ( |73| ) arises corresponds to the singularity of the reduced 
spaces of the drift system, and it gives a phase shift of 



b2 = lim f — ^ 

h^o+Jo (j\ - 7r32)^8(7r3 - ri)(r2 - 7r3)(r3 - TTg) 



= — = lim 



lim 



'"2 



1 



V8 ''^0+ {ji + r2)V?-2 - n h~*o+ Jo (ji - 7r3)^(r2 - 7r3)(r3 - tts) 

~~!= 1™ I — — 

X arctan { ^ ^^i^^^ - - + ^^^'^ - nir2 + r^)' 
\2 V(r-3 - - r2)(r2 - 7r3)(r3 - 7r3) 

-/i ji -1 -TT 



d7r3 



7r3=0 



^/8 2JlV2n^ 2 



TT 

4' 



so that the phase shift as the homoclinic orbit is approached through positive 
h is 




The same calculation for h negative gives (pi = —tt/2 and (/)2 = 7r/4 for a total 
phase shift of = 7r/2. 
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